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ABSTRACT 

A  numerical  investigation  of  unsteady  free  convection  MHD  nanofluid  flow  and  heat  transfer  over  a  stretching  vertical 
flat  plate  embedded  in  a  porous  medium  is  done.  Magnetic  field  of  constant  strength  is  applied  in  transverse  direction. 
With  the  help  of  similarity  variables,  the  governing  equations  are  sohed.  Skin-friction  coefficient  and  Nusselt  number 
are  calculated  in  non-dimensional  form.  The  effect  of  various  physical  parameters  including  in  governing  equations  are 
shown  by  graphs. 
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NOMENCLATURE 

(pCP)f 

heat  capacity  of  Uuid 

(pCP)s 

heat  capacity  of  nanoparticles 

Vf 

kinematic  viscosity  of  Uuid 

(pCp)nf 

heat  capacity  of  nanofluid 

Pnf 

NanoUuid  viscosity 

Mf 

Dynamic  viscosity  of  base  fluid 

pnf 

Density  of  nano-particles 

tm 

Free  stream  temperature 

uw 

stretching  velocity 

M 

Magnetic  number 

Gr 

Grash  of  number 

Nu 

Nusselt  number 

B0 

Applied  megnetic  field 

Pr 

Prandtl  number 

S 

Unsteady  parameter 

u 

velocity  component  along  x  direction 

K 

Porosity  of  fluid 

V 

velocity  component  along  y  direction 

T 

Temperature 

h 

similarity  variable 

4> 

Volume  fraction  of  solid  particles 

0 

Non  dimensional  tempareture 

P 

Coefficient  of  thermal  expansion 

Pf 

Density  of  Uuid 

K„f 

Thermal  conductivity  of  nanoparticles 

ps 

Density  of  solid 

Ks 

Thermal  conductivity  of  base  Uuid 

Pf 

Density  of  Uuid 

Kf 

thermal  conductivity  of  Uuid 

1.  INTRODU CTION 

Now-a-days  in  industries,  auto  mobiles  and  thermal  reactors  water,  ethylene,  oil  and  glycol  are  used  as  convectional 
Uuids.  Thermal  conductivity  of  these  lluids  is  low  which  cause  high  manufacturing  cost  and  maintenance.  In  heat 
transfer  area  traditional  fluids  are  replaced  by  nanoUuid  because  of  their  brilliant  characteristics  of 
enhanced/reduced  heat.  The  application  area  of  MHD  nanoUuid  free  convection  flow  and  heat  transfer  of  nanoUuid 
through  vertical  stretching  sheet  covers  industrial,  scientillc  and  engineering  researches.  The  new  field  of 
application  of  fluid  through  stretching  sheet  is  geothermal  energy.  Smaller  size  radiators  and  high  performance  can 
be  maintained  using  nanoUuid  as  coolants.  Heat  transfer  rate  can  be  improved  by  using  MHD  fluids.  Extrusion, 
plasma  studies,  hot  rolling,  wire  drawing,  filaments,  polymer  sheet  production,  spinning,  melt  and  production  of 
rubber  and  plastic  items  are  some  application  field  of  nanofluid  on  vertical  plate  according  to  Shankar  and  Yirga 
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(2013).  During  solution  of  viscous  fluid  over  plate  imposing  stretching  and  velocity  which  is  not  fixed,  Crane  (1970) 
introduced  the  concept  of  stretching  surface.  Gorla  and  Zinalabedini  (1987)  investigated  non-similar  free  convection  mode 
of  energy  transfer  on  vertical  plate  applying  variable  temperature  conditions.  A  numerical  solution  of  energy  process  in 
rectangular  enclosures  which  are  heated  partially  and  filled  with  nanolluid  is  solved  by  Oztop  and  Nada  (2008).  In  a  series 
of  numerical  solution  of  laminar  flow,  Khan  and  Pop  (2010)  contributed  for  flat  sheet.  Regarding  thermal  radiation  effects 
mixed  convection  flow  and  heat  transfer  over  sheet  embedded  in  porous  was  explored  by  Mukhapadhyay  (2009).  Fang  et 
al.  (2010)  examined  MHD  viscous  nanoHuid  llow  over  a  sheet  when  permeable  shirking  is  applied.  The  closed  form 
solution  of  nanoHuid  flow  past  a  stretching  sheet  is  solved  by  Hassani  et  al.  (2011)  in  context  of  boundary  layer.  Makinde 
and  Aziz  (2011)  worked  out  on  the  problem  of  nanofluid  past  a  stretching  sheet  by  applying  convective  boundary 
conditions.  By  assuming  transverse  magnetic  field,  the  solution  of  natural  convection  nanofluid  How  and  heat  transfer  has 
been  found  out  by  Hamad  (2011)  in  analytic  form.  Chamkha  et  al.  (2011)  addressed  the  effect  of  heat  generation  and 
melting  on  stretching  sheet  in  horizontal  direction  for  unsteady  nanoHuid.  Khan  and  Pop  (2011)  solving  a  similarity  form 
solution  of  nanolluid  flow  for  horizontal  plate  saturated  in  porous  medium  regarding  boundary.  In  this  series,  Ferdows  et 
al.  (2012)  worked  out  on  MHD  mixed  convective  flow  over  surface  under  exponentially  stretching.  The  effect  of  viscous 
chemical  reaction  and  dissipation  on  governing  field  of  heat  and  mass  over  stretching  surface  in  account  of  MHD  nanoAuid 
are  visualized  graphically  by  Yohannes  and  Shankar  (2013).  Vendabie  and  Sarojamma  (2014)  elaborated  the  nature  of  flow 
through  stretching  surface  under  consideration  of  magnetic  field  and  heat  generation  assuming  How  to  be  unsteady 
nanolluid.  The  numerical  solution  of  nanofluid  How  and  heat  transfer  over  exponentially  stretching  sheet  for  3-D  taking 
H2O  as  base  Huid  is  investigated  by  Nadeem  et  al.  (2014).  Yirga  and  Tesfay  (2015)  analyzed  graphically  the  effect  of 
various  parameters  including  magnetic  field  parameter  and  Soret  effects  on  nanofluid  mass  and  heat  transfer  011  permeable 
stretching  sheet.  Mahatha  et  al.  (2015)  worked  out  for  incompressible  and  viscous  over  a  stretching  surface  for  nanoHuid. 
They  have  taken  MHD  as  electrically  conducting  fluid  assuming  partially  slip  velocity,  nonlinear  thermal  radiation  and 
Newtonian  heating.  In  view  of  convective  heating  boundary  condition,  Ibrahim  and  Haq  (2016)  studied  nanoHuid  llow  of 
surface  assuming  stagnation  point. 

In  this  paper,  MF1D  nanolluid  flow  over  a  vertical  plate  in  context  of  unsteady  free  convection  embedded  in 
porous  medium  under  consideration  of  variable  temperature  is  investigated.  A  magnetic  field  of  constant  strength  is 
applied  in  transverse  direction.  The  governing  equations  are  solved  in  the  form  of  stream  function  using  similarity 
variables.  The  heat  U'ansfer  coefficient  and  skin-friction  coefficient  is  calculated  in  tabular  form.  The  effect  of  various 
parameters  included  in  governing  field  is  shown  through  graphs. 

2.  THE  PHYSICAL  DESCRIPTION  OF  PROBLEM 


The  study  for  flow  and  energy  transfer  through  a  plate  saturated  in  a  porous  medium  has  been  found  out  under 
consideration  of  stretching  and  unsteadily  in  nature.  The  x-axis  is  taken  in  vertical  direction  along  plate  and  y-axis  normal 

bx  ,  .  ,  „  ,  „  cx 


to  the  plate.  Stretching  velocity  Uv, 


1  -  at 


is  used  and  temperature  is  defined  as  T„ 


1  -  at 


.  The  plate  is  of  free  length 


vertical  direction,  so 
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Figure  1:  Physical  Diagram  of  the  Problem. 

2.1  Physical  diagram  of  the  Problem 

Applying  Bossinesq’s  approximation  and  prescribed  boundary  conditions,  the  governing  equations  are  converted  as 
The  continuity  equation 

^  =  0 
8y 

N-S  equation  of  momentum 


(1) 


du  du  1 

■  v —  = 


„  -  Mnfdprj  + 8 {pP)nf iT -Tx)~ /V 

dt  8y  pnfl  J  8y  K 


(2) 


The  governing  equation  of  temperature  field  is 


3T  dT  _  tcnf  d2T 
dt+Vdy  [pCp)nf  dy2 

The  imposed  boundary  conditions  are  given  as: 

bx 


(3) 


T  =  0 :  v  =  0,  u  =  Uw  = 


1  —  at 


(Stretching  velocity) 


r  =  r  = 


cx 


1  -at 

y  — y  oo ;  v  — y  0,  u  — ^  0,  T  — ^  T) 


(4) 


3.  MATHEMATICAL  DISCRIPTION  OF  THE  PROBLEM 

Temperature  and  Stream  function  are  used  as  given  by  Shankar  and  Yirga  (2013)  for  changing  governing  equations  in 
dimensionless  form  using  stream  function 
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y/=^vfxUwf{rj\  0(Tj)  =  ^r-^-,  rj=W^\ 

y  T«~T*  Vfx  J 

The  various  characteristics  of  nanoAuid  are  used  as  defined  by  Oztop  and  Nada(2008) 


(5) 


Pnf  =(]-<!>)Pf  +  <!>Ps 

(cpp)4  =  b-tYcpp)f+ftcpP\ 
{pP)nf={i-<fipP)f+<t>{pP\ 


»«f  = 


Pt 


(i  -<t>) 


2.5 


(Maxwell-model)(9) 


(6) 

(7) 

(8) 


(9) 


The  Maxwell-Garnetts  model  (1962)  is  used  for  the  effective  thermal  conductivity  of  nanoHuid  defined  as 


Ks  +  2 Kf  -  2(Kf  -  Ks  )(f> 
Ks  +  2k f  +  (Kf  -  Ks  )<f> 


(10) 


Inserting  the  similarity  transformation  defined  by  (5)  and  nanoHuid  physical  properties  as  given  in  equations  (6)- 
(10)  into  the  equation  (2),  energy  equation  (3)  and  corresponding  conditions  given  by  (4).  The  dimensionless  form  of 
governing  equations  are  converted  respectively  as 


r  ?  t  » _ 

7  “  ~AB 


'+(S  +  ABP  +  M 2  )f '- CAGrd 


S  —  0'—f0'  = 
2 


ED 

Pr 


0 


» t 


The  dimensionless  forms  of  boundary  conditions  are: 


(11) 


(12) 


7  =  0:  /(17)  =  0,  f\rj)  =  1,  %)  =  l| 

rj^oo:  /-> 0,  /'->0,  0^0  J 

Included  physical  parameters  are  listed  below. 


5 


—  (unsteady  parameter), 
b 


M2 


oB\ 

PfUH 


x  (magnetic  number) 


Gr  =  pf  (Tv  -  Tt  )  jt-  (Grashof  number), 

^  W 


P  = 


KU 


x  (porous  mcdium  parameter) 


(13) 
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j.  C  =  {\-<p)  +  (p- 

K,+2Kf-2<p(Kf-Ks) 

Ks+2Kf+<p(Kf-Ks) 


4.  SOLUTION  DESCRIPTION S 


(. PP)f 


Shooting  method  is  used  for  solving  governing  equations  (11)  and  (12).  In  this  method,  unknown  boundary  conditions  are 
identihed  first  and  assumed  values  are  considered.  In  this  way,  the  conditions  prescribed  on  boundary  are  transferred  to 
conditions  on  initial  value.  Then  R-K  method  is  used  to  improve  unlcnown  initial  boundaries  and  improved  continuously. 
Final  boundary  conditions  are  approached  by  solving  IVP  over  the  entire  domain. 

4.1  Solution  for  Velocity  Components 

Both  the  velocity  components  are  defined  as 


cnp_=_  I  bvf 

dx  y  (1  -  at ) 


(14) 


4.2  Heat  Transfer  Coefficient 

The  local  heat  transfer  coefficient  (Nu)  in  non-dimensional  form  is  given  as 


Nu,  = 


K,(T.,-T,y 


:.  Nu.=-E 


4.3  Skin  Friction  Coefficient 

Skin  friction  is  given  as 


C/  = 


f  du^ 
dy 


V  y  y  y= o 

The  non-dimensional  form  of  coefficient  of  slcin  friction 


(15) 


(16) 


Cf  = -2AB —j^=  f  "  (0)  (17) 

VRex 

5.  RESULTS  AND  DISCUSSIONS 

The  graphs  for  velocity  and  temperature  field  over  a  stretching  sheet  are  drawn  for  H:0-Cu  nanoHuid.  The  effect  of 
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parameters  introduced  in  momentum  and  energy  equations  are  depicted  by  graphs.  One  variable  is  changed  and  remaining 
are  fixed  and  given  specihed  values. 

5.1  Effects  on  Velocity  Profiles 

It  is  presented  by  figure  (2)  that  there  is  a  hindrance  to  fluid  flow  when  jump  in  solid  volume  fraction  of  nanoparticles 
occur  that  means  f '  (?/)  is  decreased.  Figure  (3)  portrays  the  effect  of  Grash  of  number  (Gr)  on  velocity  profile  which  is 
growing  due  to  high  temperature.  It  can  be  observed  that  behavior  of  fluid  velocity  is  reciprocal  to  the  porous  medium 
parameter  P  as  depicted  by  figure  (4).  As  the  strength  of  magnetic  force  is  increased,  it  causes  loss  in  flow  transfer  due  to 
opposite  direction  force  named  as  Lorentzian  force  and  same  is  conveyed  in  figure  (5).  The  effects  are  more  dominated  for 
powerful  magnetic  force.  The  effect  of  unsteady  parameter  on  velocity  distribution  f'(r/)  are  shown  by  figure  (6). 


Figure  3:  Effect  of  Gr  on  Yelocity  /  '(]/) 
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Figure  4:  Effect  of  P  on  Yelocity  /  '  (t]) 


11 

Figure  5:  Effect  of  M  on  Yelocity  /  '  (//) 


Figure  6:  Effect  of  S  on  Yelocity  f'(j]) 
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5.2  Effects  on  Temperature  Profiles 

The  figure  (7)  shows  enhance  heat  characteristic  of  nanolluid  which  state  that  temperature  distribution  0(r])  will  be 
increased  for  (j)  The  figure  (8)  conveys  that  higher  value  of  Gr  reduces  heat  transfer.  The  trend  of  M  and  P  is  same  for 

temperature  distribution  as  demonstrated  in  the  figure  (9)-(  10)  which  is  same  as  of  Poornima  and  Reddy  (2013).  It  is 
presented  through  the  figure  (11)  that  more  energy  is  supplied  for  higher  value  of  S.  It  can  be  seen  by  the  figure  that  heat 
transfer  is  minimum  when  unsteady  parameter  is  vanished.  The  increased  value  of  Pr  causes  less  thermal  diffusivity  due  to 
which  energy  is  increased  as  in  figure  (12). 


'n 

Figure  7:  Effect  of  (j)  on  Temperature  0(i]) 


Impact  Factor  (JCC):  8.8746 


SCOPUS  Indexed  Journal 


NAAS  Rating:  3.11 


Free  Convection  Nanojluid  Flow  Over  a  Stretching  Plate  Embedded  in 
MHD  Porous  Medium:  A  Numerical  Imestigation 


925 
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Figure  9:  Effect  of  P  on  Temperature  0(i]) . 


11 

Figure  10:  Effect  of  M  on  Temperature  0(i]) 


11 

Figure  11:  Effect  of  S  on  Temperature  0(i]) 
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11 

Figure  12:  Effect  of  Pr  on  Temperature  0(t]) 


5.3  Effects  on  Skin  Friction  Coefficient 

The  effect  of  different  physical  parameters  including  P,  <j> ,  Gr,  M  and  S  for  various  nanolluids  termed  as  AI2O3-  HiOTTi- 
HiOAg-HiO  on  skin  friction  coefficient  are  tabulated  in  Table  1.  Drag  force  increases  with  increment  in  magnetic  number 
because  of  Lorentz  drag  force.  From  table  1,  it  can  be  observed  that  HiO-Ag  and  H2O-ALO3  nanoAuid  has  highest  and 
lowest  skin  Mction  coefficient  respectively.  So  H^O-Ag  nanoHuid  has  high  skin  friction  comparatively  to  water  Aluminia 
nanolluid.  Table  1  also  presents  that  there  is  proportional  relation  between  Cf,  ((/>)  Pr,  Gr,  S  and  P.  It  is  distinguished  that 

H2O-ALO3  has  highest  skin-friction  in  respect  of  all  included  parameters. 

5.4  Effects  on  Heat  Transfer  Coefficient 

The  effect  of  different  physical  parameters  including  P,  <j> ,  Gr,  M  and  S  for  various  nanoHuids  termed  as  AI2O3-  H20,Cu- 
H^O^Ag-HiO  on  Nusselt  Number  are  tabulated  in  Table  2.  It  is  reported  from  the  table  2  that  there  is  reciprocal  relation 
between  Nussult  number  porosity  parameter  (P)  and  magnetic  number  (M).  The  increased  value  of  Pr  causes  less  thermal 
diffusivity  due  to  which  energy  is  increased  or  Nusselt  number  has  higher  value.  Nusselt  number  also  has  decreasing  trend 
in  increasing  direction  of  S  for  nanoHuids  H2O-AI2O3,  ILO-Cu  and  HzO-Ag.  The  behavior  of  Gr  and  Nu  have  same  trend. 
Due  to  high  thermal  capacity  of  silver.  HiO-Ag  has  highest  heat  transfer  rate. 

6.  CONCLUSIONS 

In  the  presence  of  porous  medium  momentum  and  energy  transfer  through  unsteadily  plate  imposing  stretching  and 
unsteadily  nature  theoretical  and  graphical  analysis  has  been  done.  The  below  mentioned  remarkable  conclusions  are 
drawn: 

•  Fluid  flow  is  of  decreasing  nature  with  increment  in  value  of  P,  M  and  (j) 

•  Fluid  flow  and  Gr  has  same  nature. 

•  Fluid  How  is  increasing  for  S  upto  tj  =0.55(approx.)  and  for  higher  value  has  reverse  effect. 
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•  Temperature  flow  and  P,  M,  Pr,  S  and  ([>  has  same  behavior. 

•  Gr  has  reverse  effect  on  heat  transfer  rate. 

•  There  is  proportional  relation  between  Drag  force  and  M. 

•  Nu  decreases  with  increment  in  amount  of  M,  P,  Pr,  S  and  (f>  . 
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APPENDIEX 


Table  1:  Effect  of  Physical  Parameters  on  Coefficient  of 
Skin-Friction  (Cf)  of  Different  Type  of  Nanolluids 


Physical  Parameters 

Nanolluid 

H2O-AI2O3 

H2O-CU 

H20-Ag 

M 

Pr 

P 

<l> 

Gr 

S 

Cf 

Cf 

Cf 

0 

5.784 

1 

0.04 

2 

1.5 

1.253158 

1.32489 

1.345468 

0.5 

5.784 

1 

0.04 

2 

1.5 

1.332998 

1.416227 

1.440311 

1 

5.784 

1 

0.04 

2 

1.5 

1.560647 

1.675028 

1.708553 

2 

5.784 

1 

0.04 

2 

1.5 

2.335491 

2.541652 

2.602642 

0.5 

0.7 

1 

0.04 

2 

1.5 

1.343393 

1.42776 

1.452273 

0.5 

1 

1 

0.04 

2 

1.5 

1.342828 

1.427091 

1.451623 

0.5 

1.5 

1 

0.04 

2 

1.5 

1.341817 

1.425973 

1.450458 

0.5 

2 

1 

0.04 

2 

1.5 

1.340801 

1.42485 

1.449289 

0.5 

5.784 

1 

0.04 

2 

1.5 

1.332998 

1.416227 

1.440311 

0.5 

5.784 

0 

0.04 

2 

1.5 

1.003724 

1.09775 

1.125056 

0.5 

5.784 

1 

0.04 

2 

1.5 

1.332998 

1.416227 

1.440311 

0.5 

5.784 

2 

0.04 

2 

1.5 

1.62991 

1.704693 

1.726253 

0.5 

5.784 

3 

0.04 

2 

1.5 

1.900896 

1.968923 

1.988464 

0.5 

5.784 

1 

0.04 

2 

1.5 

1.332998 

1.416227 

1.440311 

0.5 

5.784 

1 

0.1 

2 

1.5 

1.574072 

1.790115 

1.850355 

0.5 

5.784 

1 

0.04 

1 

1.5 

1.628536 

1.713603 

1.739026 

0.5 

5.784 

1 

0.04 

2 

1.5 

1.332998 

1.416227 

1.440311 

0.5 

5.784 

1 

0.04 

3 

1.5 

1.043372 

1.124684 

1.147499 

0.5 

5.784 

1 

0.04 

4 

1.5 

0.759304 

0.838629 

0.860241 

0.5 

5.784 

1 

0.04 

2 

0 

1.069526 

1.094935 

1.102844 

0.5 

5.784 

1 

0.04 

2 

1 

1.248805 

1.313378 

1.33232 

0.5 

5.784 

1 

0.04 

2 

1.5 

1.332998 

1.416227 

1.440311 

0.5 

5.784 

1 

0.04 

2 

3 

1.589226 

1.724743 

1.763609 

0.5 

5.784 

1 

0.04 

2 

5 

1.974349 

2.166915 

2.223248 

Table  2:  Effect  of  Physical  Parameters  on  Heat  Transfer  Coefficient 
(Nusselt  Number)  Nux  of  Different  Type  of  Nanoduitls 


Physical  Parameters 

Nanolluid 

H2O-AI2O3 

H2O-CU 

H20-Ag 

M 

Pr 

P 

<!> 

Gr 

S 

Nu 

Nu 

Nu 

0 

5.784 

1 

0.04 

2 

1.5 

0.973095 

0.966409 

0.964098 

0.5 

5.784 

1 

0.04 

2 

1.5 

0.965023 

0.957399 

0.954696 

1 

5.784 

1 

0.04 

2 

1.5 

0.942609 

0.932643 

0.928947 

2 

5.784 

1 

0.04 

2 

1.5 

0.873271 

0.85837 

0.852401 

0.5 

0.7 

1 

0.04 

2 

1.5 

0.995795 

0.994997 

0.994522 

0.5 

1 

1 

0.04 

2 

1.5 

0.994148 

0.992839 

0.992378 
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0.5 

1.5 

1 

0.04 

2 

1.5 

0.991189 

0.989225 

0.988532 

0.5 

2 

1 

0.04 

2 

1.5 

0.988209 

0.985588 

0.984663 

0.5 

5.784 

1 

0.04 

2 

1.5 

0.965023 

0.957399 

0.954696 

0.5 

5.784 

1 

0.04 

2 

1.5 

0.965023 

0.957399 

0.954696 

0.5 

5.784 

2 

0.04 

2 

1.5 

0.935969 

0.92988 

0.927293 

0.5 

5.784 

3 

0.04 

2 

1.5 

0.910819 

0.905952 

0.903435 

0.5 

5.784 

1 

0 

2 

1.5 

0.975324 

0.975324 

0.975324 

0.5 

5.784 

1 

0.04 

2 

1.5 

0.965023 

0.957399 

0.954696 

0.5 

5.784 

1 

0.1 

2 

1.5 

0.950415 

0.935306 

0.929274 

0.5 

5.784 

1 

0.04 

1 

1.5 

0.933962 

0.926684 

0.923621 

0.5 

5.784 

1 

0.04 

2 

1.5 

0.965023 

0.957399 

0.954696 

0.5 

5.784 

1 

0.04 

3 

1.5 

0.995133 

0.987186 

0.984825 

0.5 

5.784 

1 

0.04 

4 

1.5 

1.024324 

1.016074 

1.014037 

0.5 

5.784 

1 

0.04 

2 

0 

1.667883 

1.658433 

1.66351 

0.5 

5.784 

1 

0.04 

2 

1 

1.198691 

1.190114 

1.189972 

0.5 

5.784 

1 

0.04 

2 

1.5 

0.965023 

0.957399 

0.954696 

0.5 

5.784 

1 

0.04 

2 

3 

0.387218 

0.383754 

0.377315 

0.5 

5.784 

1 

0.04 

2 

5 

0.063236 

0.06302 

0.060218 
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